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Three distinct scales are identified in the excitation spectrum of the gluon field 
around a static quark-antiquark pair as the color source separation R is varied. 
The spectrum, with string-like excitations on the largest length scales of 2-3 fm, 
provides clues in its rich fine structure for developing an effective bosonic string 
description. New results are reported from the three-dimensional Z(2) and SU(2) 
gauge models, providing further insight into the mechanism of bosonic string for- 
mation. The precocious onset of string— like behavior in the Casimir energy of the 
static quark-antiquark ground state is observed below R=l fm where most of the 
string eigenmodes do not exist and the few stable excitations above the ground 
state are displaced. We find no firm theoretical foundation for the widely held 
view of discovering string formation from high precision ground state properties 
below the 1 fm scale. 



'Speaker at the conference. 
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1. QCD String Spectrum and the Casimir Energy 

Last year, we presented a new analysis of the fine structure in the QCD 
string spectrum at the Lattice 2002 conference. Shortly afterwards, two 
papers were submitted using complementary methods for finding definitive 
signals of bosonic string formation from the rich excitation spectrum^ and 
the ground state Casimir energy^ 

QCD String Spectrum 

Three exact quantum numbers which are based on the symmetries of 
the problem determine the classification scheme of the gluon excitation 
spectrum in the presence of a static qq pair.^ We adopt the standard no- 
tation from the physics of diatomic molecules and use A to denote the 
magnitude of the eigenvalue of the projection Jg • R of the total angular 
momentum Jg of the gluon field onto the molecular axis with unit vector 
R. The capital Greek letters S, 11, A, . . . are used to indicate states with 
A = 0, 1, 2, 3. ... , respectively The combined operations of charge conju- 
gation and spatial inversion about the midpoint between the quark and the 
antiquark is also a symmetry and its eigenvalue is denoted hy ijcp- States 
with r]cp = 1(— 1) are denoted by the subscripts g{u). There is an addi- 
tional label for the S states; E states which are even (odd) imder a reflection 
in a plane containing the molecular axis are denoted by a superscript + (— )• 
Hence, the low-lying levels are labeled E+, E~, S+, E~, Ilg, n„, Ag, A„, 
and so on, E+ designating the ground state. For convenience, we use T to 
denote these labels in general. For better resolution of the fine structure 
in the spectrum, the gluon excitation energies Er(R) were extracted from 
Monte Carlo estimates of generalized large Wilson loops on lattices with 
small at/as aspect ratios and improved action. Restricted to the R=0.2-3 
fm range of a selected simulation, the energy spectrum is shown in Fig. 1 
for 10 excited states. 

On the shortest length scale, the excitations are consistent with short 
distance physics without string like level ordering in the spectriim. A 
crossover region below 2 fm is identified with a dramatic rearrangement 
of the short distance level ordering. On the largest length scale of 2-3 fm, 
the spectrum exhibits string- like excitations with asymptotic tt/R string 
gaps which are split and slightly distorted by a fine structure. It is remark- 
able that the torelon spectrum of a closed string, with one unit of winding 
number around a compactified direction, exhibits a similar fine structure 
on the 2 3 fm scale, as reported for the first time at Lattice 2003.'^ This 
finding eliminates the boundary effects of fixed color charges as the main 
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Figure 1. Short distance degeneracies which are not string-like and their crossover to- 
wards the QCD string spectrum are shown from Ref. 1 where further details are ex- 
plained. The color coded solid curves with simulation points, which identify energy 
levels degenerate in the asymptotic string limit, are only shown for visualization and do 
not represent fits to the data. The yellow line without data points marks a lower bound 
for the onset of mixing effects with glueball states which requires careful interpretation. 
The symbol LW indicates the R range of high precision Casimir energy calculations from 
Ref. 2. 



source of the fine structure in the distorted spectrum. 
Casimir Energy 

In a complementary study,^ a string-hke Casimir energy and the re- 
lated effective conformal charge, Cofr(R) = — I2R^F'(R)/(7r(D — 2)), were 
isolated where F(R) is the force between the static color sources and D 
is the space-time dimension of the gauge theory with bosonic string for- 
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mation. With unparalleled accuracy, Ceff(R) was determined for the gauge 
group SU(3) in D=3,4 dimensions, in the range 0.2 f m < R < 1.0 fm, below 
the crossover region of the string spectrum. A sudden change with increas- 
ing R, well below 1 £m, was observed in Ceff(R), breaking away from the 
the short distance running Coulomb law towards the string- like Ccff (R) ~ 1 
behavior. This was interpreted as a signal for early bosonic string forma- 
tion. The results are surprising because the scale R is not large compared 
with the expected width of the confining flux, and more quantitatively, the 
string-like Casimir energy behavior is observed in the R range where the 
spectrum exhibits complex non-string level ordering, as shown in Fig. 1. We 
will try to develop now a better understanding of the seemingly paradoxical 
situation. 

2. Bosonic String Formation in the Z(2) Gauge Model 

The three-dimensional Z(2) lattice gauge model represents considerable 
simplification in comparison with four-dimensional lattice QCD. The SU(3) 
group elements on links are replaced by Z(2) variables, and the reduction 
to three dimensions implies a nontrivial continuum limit with a finite fixed 
point gauge coupling. 

Dual 4>^ Field Theory Representation 

The main features are easily seen from the dual transformation of the 
Z(2) gauge model to Ising variables which can be replaced by the real scalar 
field of 0* field theory in the critical region.*^ The continuum model exhibits 
confinement and bosonic string formation in the broken phase of the Ising 
representation. In addition, a nontrivial glueball spectrum is observed^ 
with finite masses when measured in units of the string tension. The string 
tension of the confining flux in the Z(2) gauge model becomes the inter- 
face energy in the dual Ising-(/)^ representation, and the lowest mass 0"*" 
glueball state of the gauge theory with mass m maps into the elementary 
scalar of the dual lattice, with inverse correlation length m in the criti- 
cal region. Higher glueball states are Bethe-Salpeter bound states of the 
elementary scalar.^ The dual (h^ Lagrangian of the Z(2) gauge model, in 
rigorous theoretical setting, is in analogy with the dual Landau-Ginzburg 
superconductivity model which attempts to describe the unknown micro- 



^Rcfcrcnccs to earlier work on the three-dimensional Z(2) gauge model can be found in 
a recent paper on the finite temperature properties of the Z(2) string.* 
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scopic quark confinement mechanism of QCD. ^ 

The Ising-</)^ field theory model is particularly intriguing from the micro- 
scopic string theory viewpoint, if we recall Polyakov's work on the connec- 
tion with the theory of random surfaces. Using loop equations of closed Wil- 
son loops near the continuum limit, he conjectured the equivalence of the 
three-dimensional Z(2) lattice gauge theory to a fermionic string theory/ 

Renormalization Scheme 

In three euclidean dimensions, the Z(2) model is described in the crit- 
ical region (continuum limit) by a real order parameter field (j) with the 
Lagrangian 

c=-kd,M'-^i<^i-^-^r. (1) 

^ 4! go 

The most frequently used renormalization scheme requires in the broken 
phase that the tadpole diagrams completely cancel without coupling con- 
stant renormalization {g = go) and with wave function renormalization 
(t>a = \/Z(l). In the following, with lattice cutoff, we define a scheme with 
finite coupling constant renormalization keeping the renormalized mass of 
the elementary scalar exactly at the pole of its propagator. Since the wave 
function renormalization is finite to every order, for convenience we choose 
Z—1 in 1-loop calculations. With go = g -I- Sg, Vq = Smp/go = v^ + (5v^, 
the renormalized Lagrangian for elementary excitations rj around the vac- 
uum expectation value (p = v is the starting point of the renormalized loop 
expansion with two counterterms to one-loop order, 

5v2=limv^oo ^g-^E^o > ^^=^'^- (2) 

k k 

The infinite spatial volume limit is taken in the sum over the spectrum of in- 
verse lattice energies of free massive excitations with periodic boundary 
conditions. The coupling constant counterterm Sg satisfies the renormal- 
ization condition on the physical propagator pole to one-loop order. 

In the presence of a pair of static sources, represented by a Wilson 
loop in the Z(2) gauge model, the renormalization procedure is unchanged 
in the Ising-^'^ field theory description. The only change is in the lattice 
Lagrangian where the sign of the nearest neighbor interaction term is flipped 
on links which piuicture the surface of the Wilson loop on the dual Z(2) 



''A recent review of quark confinement and dual superconductivity is given in Ref. 6 
with discussion of earlier work and references. 
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gauge lattice. This flip represents a disorder line, or seam, between the two 
static sources on the spatial lattice. The end points of the seam are fixed 
but otherwise it is deformable by a "gauge transformation" of variables 
without changing the partition function. This invariance is inherited from 
the gauge invariant representation of the Wilson surface in the Z(2) gauge 
model. 

Numerical Implementation of the Loop Expansion 

The dual transformation of the Z(2) model to Ising variables facilitates 
very efiicient simulations with multispin coding. The loop expansion pro- 
vides theoretical insight into Monte Carlo simulations of the excitation 
spectrum using high statistics multispin Ising codes complemented by cj)'^ 
field theory codes. Since the fixed point value u* of the dimensionless cou- 
pling constant u = g/m is not small, the simulations provide an important 
cross-check on the convergence of the loop expansion which itself has to be 
implemented in a numerical procedure. The renormalized loop expansion 
in the presence of static sources requires the following three steps. 

(i) First, for a given physical mass m and renormalized coupling g, the time 
independent renormalized classical field equation, 

a? -m <^,(x,y) + -</.,(x,y) =0 , (3) 

of the static soliton (f)s is solved on the lattice in the f* = (x, y) plane with 
flipped nearest neighbor interaction links along the seam between static 
sources. In the Z(2) gauge model, the two sources can be interpreted as 
opposite sign charges with an electric flux connecting them. In the (f>'^ 
representation we refer to the (ps classical solution as a static soliton, rather 
than the earthy flux-tube term. In the numerical procedure, a generalized 
Newton type nonlinear iterative scheme was implemented to obtain to 
double precision accuracy. 

(ii) Second, the fluctuation spectrum around the static soliton 0s(x, y) is 
determined by splitting the field into the classical solution plus fluctua- 
tions, (/)(x, y,t) = 0s(x, y) -|- 77(x, y,t), with the eigenmodes of the fluctua- 
tion field r/(x,y,t) = En[an(t)V'n(x, y) + al{t)ip*{x,y)] satisfying the eigen- 
value equation 

The time dependence of the fluctuation fleld rj is given in interaction pic- 
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ture by an(t) = an(0)e~'^"* where the Hamihonian is spUt into a quadratic 
part and an interaction part of the ipn eigenniodes. The second derivative 
of the IJ{4>) = j!{4>'^ — ^^)^ renormahzed field potential energy is taken 
with respect to (f) in Eq. (4), with cj) = (j)g substituted subsequently. Two 
parity quantum numbers Px- Pj- split the cigcnmodcs into four separate 
symmetry classes. With the two sources located at (x, y) = (R/2,0) and 
(x, y) = (— R/2, 0), the quantum number Fx = ±1 corresponds to the reflec- 
tion symmetry x — > — x of V'n(x, y) and Py = ±1 corresponds to the y — > — y 
reflection symmetry. The full spectrum of eigenvalues and eigenfunctions 
of Eq. (4) are computed by an Arnold! diagonalization procedure in the 
flnite volume of the lattice. Using the parity symmetries of the theory, 
diagonalization of large lattices with sizes up to 200x200 in the (x,y) plane 
were performed. 

(iii) Third, the systematic renormahzed loop expansion with the (j)a static 
soliton background is developed by building the finite volume field propa- 
gator DF(r, t;r',t') in Minkowski time, 

n ^ 

from the static Vn(f) eigenmodes. An euclidean rotation is performed on 
the propagator during the numerical evaluation of the loop diagrams. The 
counterterms 5v^ and 5g are used to remove loop divergences in the contin- 
uum limit and to keep the exact pole location at the physical mass m. Using 
the propagator of Eq. (5), the fluctuation correction to the static soliton 
profile (/>s was calculated to one-loop order, together with similar calcula- 
tions of the ground state energy and excitation energies. In this work, we 
only report numerical results on the fluctuation spectrum of Eq. (4) and its 
1-loop contribution to the ground state energy. 

String Excitations in the Loop Expansion 

For sufficiently large R, the discrete Py = — 1 bound state spectrum of 
Eq. (4) is expected to evolve into the asymptotic En = ttN/R (N = 1,2,.. .) 
Dirichlet string spectrum of massless string excitations which originate from 
the translational mode of the well-known one-dimensional <j)^ soliton by the 
following simple consideration. 

Consider first the spatial lattice in the finite (x,y) plane with a seam of 
flipped links winding around the compact x-direction with periodic bound- 
ary condition. The classical solution </>tor(y) = ui-\/3/g tanh(m|y|/2) of 
Eq. (3) defines the torelon which is independent of x and winds around the 
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compact X direction with a scam positioned at y=0. We use contimmm 
notation for the torelon and its excitations, but finite cutofi^ and volume 
effects are included in the numerical work. For x > 0, the transverse profile 
of the torelon is identical to that of the well-known one- dimensional soli- 
ton, and for x < a sign flip is involved because of the seam at y=0. The 
torelon eigenmodes of Eq.(4) with Py = — 1 have the simple form 



with quantized momenta, p^ = 27rN/R, N = ±1,±2, . . ., running along x 
in the compact interval R with periodic boundary condition. The energy 
spectrum is given by En = 27rN/R, with positive N values. 

The classical transverse profile (/>tor(y) of the torelon coincides, to a good 
approximation, with that of the static soliton if the separation between 
the sources is large enough. The static soliton profile (f)s does not interpolate 
from -v at large negative y to +v for large positive y at fixed x because of 
the flipped links along the seam. Rather, i/ig approaches v everywhere, far 
away from the seam line. The eigenmodes of the fluctuation operator are 
restricted now between the two sources and they are close to the form 



V-r^'Cx, y) = J A 'PUy) ■ VV^ sin(7rN/x), -R/2^x^ R/2 , (7) 



with N taking positive integer values. The spectrum of these standing waves 
is the same as that of a massless Dirichlet string oscillating in the (x,y) 
plane with fixed ends. The excited eigenmodes of the effective Schrodinger 
equation, like the one of Fig. 2b for N—S, are therefore in one-to-one cor- 
respondence with massless Dirichlet string oscillations. The spectrum and 
the wave functions are expected to be somewhat distorted at finite R be- 
cause of the distortions of the effective Schrodinger potential around the 
sources. 

Representative examples of the numerical work are shown in Fig. 2 
where the static soliton solution (/)s(x, y) on a 160x80 spatial lattice in 
the (x,y) plane corresponds to source separation R=100 and physical mass 
m=0.319 in the critical region (all dimensional quantities are expressed in 
lattice spacing units). For later comparisons, the lattice correlation length 
m~^ and the renormalized coupling g were chosen in the critical region to 
match one of our Monte Carlo simulations with v=0.45 and string tension 
a = 0.0101. The static soliton solution (/),s (f') determines the attractive po- 
tential energy of the effective Schrodinger eigenvalue problem in Eq. (4) 
which has a discrete bound state spectrum and a nearly continuous dense 
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Figure 2. The static soliton solution (f)s of field equation (3) is shown in (a) witii the 
choice R=100 and renormalized parameters given in the text; (b) shows the N=8 massless 
string excitation of the static soliton from the numerical diagonalization of the eigenvaule 
equation (4); the second massive string excitation (K=2 breathing mode) is shown in 
(c) from the numerical solution of the eigenvalue equation, and compared in (d) to exact 
Monte Carlo simulation of the same state with remarkable agreement. 

spectrum above the glueball threshold m representing scattering states on 
the static sohton in the infinite lattice volume limit. 

The one-dimcnsional soliton, with classical mass 2m*^/g, has a 
massive intrinsic excitation, or breathing mode, whose excitation en- 
ergy is (\/3/2)m. In the large R limit, the intrinsic excitations 
of (j)s with Py = +1 become massive breathing modes of the Dirich- 
let string. The asymptotic spectrum of a massive string, given by 
Ek = \/3m2/4 + 7r^K2/R2, K = 1, 2, . . ., is associated with eigenmodes 
like the K=2 wave function of Fig. 2c. The corresponding standing wave 
solutions, 

V'K(x,y) = V3m/4 sech(my/2) tanh(m|y|/2) • ^2/^ sin(7rK/x) , (8) 

originate from the massive excitations of the torelon with restriction to 
standing waves in the — R/2 ^ x ^ R/2 interval. 
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3. Probing the String Theory Limit 

We describe the t) fluctuations around the static sohton (ps by the sum 
of three fields, 77(1', t) = ^(r, t) + x{r, t) + </3(r, t), with 

^(x,y,t) = ^[aN(t)V'N(x,y) + aJj(t)V'N(x,y)] , 

N 

X(x,y,t) = ^[aK(t)V'K(x,y) + aJ^(t)V'K(x,y)] , 

K 

'y^(x,y,t) = ^[an(t)V'n(x,y) + a^(t)V'*(x,y)] , 

n 

where ^ is restricted to bound states with negative Py parity which are 
expected to evolve into massless string excitations for large R. The field x is 
restricted to Py = +1 parity bound states which evolve into massive string 
excitations, and ip \s a sum over scattering states above the 0+ glueball 
threshold m in the continuum. 

These fields are coupled in the interaction Lagrangian, and when the 
massive fields x ^-nd Lp are integrated out, we get a nonlocal Lagrangian in 
the ^(x, y, t) field describing massless string excitations in the large R limit. 
As indicated by Eq. (7), the y-dependence in all the Py = — 1 parity bound 
state wave functions is approximately factored out in the large R limit. 
Hence, the ^(x, y, t) field can be replaced on large length scales by the field 
f(x, t) which becomes the geometric string variable of low energy excitations 
measuring the displacements of the flux center-line in the y-direction as a 
function of x and t. 

Effective String Action 

The nonrcnormalizable effective action of the f(x, t) field, with the mas- 
sive fields integrated out, is given in a derivative expansion by 

s«(t) = - / dxd. (1 - \imf - \mff - + ■ ■ ■ ) 

+ const ■ {difU{dif + ■ , (9) 

where d ^ is a long wavelength expansion. The string tension cr can 
be expressed as a = {2-Ka')~^ , and the notation (3f)^ = d^id^f is used in 
Eq. (9). Since / is related to massless Goldstone excitations, originating 
from the restoration of translation invariance in torelon quantization, only 
derivatives of f appear in Scff. The first three terms in the derivative ex- 
pansion come from the kinetic terms in the original field theory action. 
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They are independent of the details of the field potential except for the 
overall factor of (27ra')~^. The first line in Eq. (9) agrees with the equiv- 
alent terms of the Nambu-Goto (NG) action, Sng = (27ra')~^-\/l — (9f)^, 
when expanded in (9f)^. The second line in Eq. (9) has contributions from 
the geometric curvature, but new contributions also appear whose geomet- 
ric origin remains unclear. This is where the effective string action begins 
to show deviations from the NG string. To correct for end effects around 
the static sources, the effective action of Eq. (9) has to be augmented by 
boundary operators for the complete description of the Dirichlet string. 

Exact Excitation Spectrum, from Simulations 

First, the physical mass, the vacuum expectation value v, and the renor- 
malized coupling g were determined in high precision Monte Carlo simula- 
tions of the bulk euclidean lattice action in the Ising-t^^ field theory rep- 
resentation without the seam of flipped links. The rcnormalizcd physical 
parameters were used as input to compare the loop expansion with simu- 
lations in the presence of static sources. The excitation spectrum around 
the static sources was determined from Monte Carlo estimates of correla- 
tion matrices which included an extended set of optimized operators. These 
operators were built from eigenfunctions of the fluctuation operator on two- 
dimensional time slices of the lattice. 

Simulation results arc shown in Fig. 3 for string formation as R is 
stretched. The shape of the static soliton profile (t>s{y^, y) around the sources 
is depicted in the (x,y) plane of the two-dimensional spatial lattice for R=10 
and R=100. With the scale set by the string tension, the two R values cor- 
respond to 0.3 fm and 3 fm, respectively. The rcnormalizcd bulk physical 
parameters of the simulations were given in the discussion of Fig. 2 as 
m=0.319, v=0.45, and a = 0.0101. At the smaller R value, the static soli- 
ton with bag like shape is not stretched, and there is only one bound state 
excitation below the glueball threshold. At R=3 fm, the stretched static 
soliton supports many string excitations, with the N=l,4,8 wave functions 
displayed. The Monte Carlo simulations are in good agreement with results 
from the eigenmodcs of the fluctuation operator. This is indicated by the 
good match of Fig. 2a versus Fig. 3c and Fig. 2b versus Fig. 3f. 

Matching the Excitation Spectrum, to the String Action 

A large number of excitation spectra were obtained from highly accurate 
simulations as R was varied in a wide range from 0.3 fm to approximately 
10 fm. A good test for string formation is provided by the behavior of the 



February 1, 2008 4:27 WSPC/Trim Size: 9in x 6in for Proceedings 



riken jk 



12 






Figure 3. The soliton profile for R=10 is shown in (a), with the only bound state wave 
function below the glueball threshold depicted in (b). At R=100, the stretched soliton 
configuration shown in (c) exhibits several string-like excitation with the N=l wave 

function shown in (d), N=4 in (c), and N=8 in (f). The simulation results match the 
loop expansion of Fig. 2 with common renormalized parameters. 

spectrum as a function of R. The NG spectrum, 

/ D- 2 27rN\^ , , 

with fixed end boundary conditions in D dimensions, was calculated in 
Ref. 8. N=0 corresponds to the string ground state and positive integer N 
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Figure 4. The energy gap AE above the ground state is plotted as AE/(N7r/R) — 1 
to show percentage deviations from the asymptotic N = 1 string level. Several Z(2) 
simulations with cyan, blue, red, and green points are combined with good scaling prop- 
erties. The open circles represent D=3 SU(2) results after readjusting the ratio of the 
string tension a to the glueball mass in Z(2). The null line coresponds to the tree level 
7r/N NG string gaps. The dashed blue and green lines are 1-loop and two-loop NG 
approximations, repectively, and the black line is the full NG prediction. 



values label the excitations of the Dirichlet string. Although there exists an 
inconsistency in the quantization of angular momentum rotations around 
the qq-axis at finite R values unless D = 26, the problem asymptotically 
disappears in the R — s- cx) limit. * This is expected from the earlier discussion 
on string formation in the loop expansion. Indeed, derived in a consistent 
D=3 0** field theory, the first few terms of the effective string action match 
the coefficients of their NG counterparts as seen in Eq. (9). If a string limit 
is reached for a large enough R range, the expansion of En into inverse 
powers of R from Eq. (9) should agree with the simulations at least to one 
nontrivial order. At small R values the expansion will break down since the 
mathematical NG string will be the inconsistent description of the bag-like 
soliton and its excitations in three dimensions. 

The comparison of simulations to string theory is illustrated in Fig. 4 
where the exact N=l excitation is plotted against the numerical spectrum 
of the fluctuation operator and the predictions of the Nambu-Goto string 
model. The numerical spectrum from Eq. (4) is the rcnormalizcd tree level 
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starting point of the loop expansion in 0"* field theory setting without as- 
suming string formation. It is close in shape and details to the exact results 
for the entire R range, as shown by the solid red line with the soliton tag 
in Fig. 4. It is expected that higher loop corrections will bring the agree- 
ment even closer. The simulations, however, deviate substantially from the 
predictions of the loop expansion in the NG string model, particularly at 
smaller R values below 1 fm where the loop expansion suddenly begins to 
diverge. The string formation for large R is clearly seen in the spectrum 
of Fig. 4 and all the other spectra we obtained, but further work is needed 
for quantitative matching of the coefficients in the effective string action. 

Our results differ from the findings of Ref. 4 where simulations of the 
finite temperature Z(2) string were reported in good agreement with the 
expansion of the NG model to 1-loop order below the R=l fm scale. 



4. The Casimir Energy Puzzle 

The breakdown of the effective string description below 1 fm and the related 
Casimir energy paradox are illustrated with the calculation of the ground 
state energy from the renormalized fluctuation operator around the static 
4>s soliton. In the renormalized loop expansion, the early onset of the 
string- like Ceff(R) « 1 behavior in the R range below 1 fm is not associated 
with massless string eigenmodes which are mostly missing for R < 1 fm. 
The results are consistent with the direct simulations of Cofr(R) and the 
spectrum in our Z(2) model. The simulations of Refs. 1, 2 present the 
same puzzle in QCD. If there is any physics associated with this puzzle, it 
remains unresolved. 

Casimir Energy from the Fluctuation Operator 
The soliton ground state energy Eg in Eq. (11), 



- A^v2 
12 





r. 3m2l 




r, 3m21 


/ 


4>s 


dxdy+l/ 


0s 




g . 




g . 



dxdy , (11) 



includes the renormalized classical energy E^', the sum of zero point energies 
summed over all eigenmodes a of the fluctuation operator of Eq. (4) , and a 
sum over all momenta k in the zero point energy of the bulk vacuum without 
sources. The difference of the two cigcnmodc sums is still divergent in 
the continuum limit. The counterterm removes this divergence and 
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contributes a finite r-dependent energy density to the ground state. The Sg 
counterterm also contributes a finite and r-dependent ground state energy 
density. 

Exact Ground State Results 

The simple definition CcS (R) = -24R2(E^(R) - cr)/(7r(D - 2)) was 
used to isolate the effective Casimir energy term in the ground state of 
the soliton. The first derivative Eg was directly simulated on the lattice by 
a special method we developed. The string tension a was determined in 
high precision separate runs from the ground state of long torelons. The 
simulations in Fig. 5 are also compared with the predictions of the Nambu- 
Goto (NG) string model and the predictions of Eq. (11) from the numerical 
evaluation of the fiuctuation operator spectrum. 

The agreement of Ceff(R) with exact simulation results, as determined 
from the fluctuation operator spectrum, is quite good down to very small 
R values. At smaller R values, full agreement can only be expected from 




0.7 I ^ ' ' ' ' ' ' ' 1 

5 10 15 20 25 30 35 40 45 



Figure 5. The red points are from high precision Ising-fli'* simulations. The solid black 
curve with NG label is the full NG prediction, Ccff (R) = 1 is the asymptotic string 
result which corresponds to the tree-level NG prediction. The dashed green line shows 
the 1-loop approximation which includes the first correction to tree level from the 
expansion of Eq. (10). The solid red line is calculated from the numerical evaluation of 
Eq. (11). The dashed blue line is obtained from the full red line by subtracting the E^' 
contribution. 



February 1, 2008 4:27 WSPC/Trim Size: 9in x 6in for Proceedings 



riken jk 



16 



-0.001 - 



0.003 r 



0.002 - 



0.001 




80 



Figure 6. The ground state fluctuation energy density is plotted for R=40. Tiie two 
peaks represent large fluctuation contributions around the static sources. 

higher loop corrections. It is interesting that the classical contribution of (f)^ 
to Ccff (R) is significant below 1 fm. At large R, in the true string formation 
limit, the higher loop corrections should not contribute to the asymptotic 
1/R Casimir term. 

The dominant contribution to the ground state energy of the static soli- 
ton from the fluctuation operator is coming from the continuum spectrum 
but the ground state energy density remains concentrated around the static 
soliton. This is also seen in the Casimir energy density which is calculated 
from the first derivative of the energy density with respect to R, integrat- 
ing to Ceff(R) as we explicitly checked. For illustration, the Casimir energy 
density is shown in Fig. 6. 

Although the 1/R expansion of Ceff(R) from the NG prediction of 
Eq. (10) is not divergent below 1 fm, it breaks away from the data in a 
rather dramatic fashion. 

5. Conclusions 

We established bosonic string formation in a large class of gauge theory 
models from a direct study of the excitation spectrum at large separation 
of the static sources. The spectrum, with string-like excitations on the 
length scales of 2-3 fm and beyond, provides clues in its rich fine structure 
for developing an effective bosonic string description. The matching of the 
string-like spectrum to an effective string action remains a challenge. 



February 1, 2008 4:27 WSPC/Trim Size: 9in x 6in for Proceedings 



riken jk 



17 

Our results at small R differ from, the findings of Ref. 4 where simulations 
of the finite temperature Z(2) string were reported in good agreement with 
the expansion of the NG model to 1-loop order below the R=l fm scale. 
This agreement was interpreted as further support at finite temperature for 
the precocious onset of bosonic string formation in QCD below the 1 fm 
scale as reported in Ref. 2. 

We find no firm theoretical foundation for discovering string formation 
from high precision ground state properties below the 1 fm scale. The ex- 
planation for string like finite temperature free energy behavior below 1 fm 
also remains unclear. Further work is needed to understand the interpreta- 
tion of the results from Refs. 2, 4 which present impressive high precision 
simulations. 
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